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We use a driven Monte Carlo dynamics in the phase representation to determine the linear re- 
sistivity and current-voltage scaling of a two-dimensional Josephson-junction array at an irrational 
flux quantum per plaquette. The results are consistent with a phase-coherence transition scenario 
where the critical temperature vanishes. The linear resistivity is nonzero at any finite tempera- 
tures but nonlinear behavior sets in at a temperature-dependent crossover current determined by 
the thermal critical exponent. From a dynamic scaling analysis we determine this critical exponent 
and the thermally activated behavior of the linear resistivity. The results are in agreement with 
earlier calculations using the resistively shunted-j unction model for the dynamics of the array. The 
linear resistivity behavior is consistent with some experimental results on arrays of superconducting 
grains but not on wire networks, which we argue have been obtained in a current regime above the 
crossover current. 
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Most theoretical investigations of the vortex-glass 
phase in superconductors have considered model systems 
where there is a combined effect of quenched disorder and 
frustration 1 . However, in artificial Josephson-junction 
arrays, frustration without disorder can in principle be 
introduced by applying an external magnetic field on 
a perfect periodic array of weakly coupled supercon- 
ducting grains^ 3 -^ and similarly on superconducting wire 
networks^. The frustration parameter /, the number of 
flux quantum per plaquette, is given by / = 4>/4> , the 
ratio of the magnetic flux through a plaquette <p to the su- 
perconducting flux quantum 4> = hc/2e. ft can be tuned 
by varying the strength of the external field. Frustration 
effects can be viewed as resulting from a competition be- 
tween the underlying periodic pinning potential of the 
array and the natural periodicity of the vortex lattice^. 
At a rational value of /, the ground state is a commensu- 
rate pinned vortex lattice leading to discrete symmetries 
in addition to the continuous U(l) symmetry of the su- 
perconducting order parameter. The resistive transition 
is only reasonably well understood for simple rational 
values of /. 

At irrational values of /, the resistive behavior is much 
less understood since the vortex lattice is now incom- 
mensurate with the periodic array, fn early Monte Carlo 
(MC) simulations^ the ground state was found to con- 
sist of a disordered vortex pattern lacking long range or- 
der which could be regarded as a some sort of vortex- 
glass state without quenched disorder. Glassy-like be- 
havior was indeed observed in these simulations suggest- 
ing a possible superconducting (vortex-glass) transition 
at finite temperatures. However, some arguments also 
suggested that the critical temperature should vanish 7 -^. 
Simulations of the current-voltage scaling using the resis- 
tively shunted-junction model for the dynamics of the ar- 
ray found that the behavior was consistent with an equi- 
librium resistive transition where the critical temperature 



vanishes^, similar to the resistive transition described by 
the the gauge-glass model in two dimension*^, but with 
different values for the correlation-length critical expo- 
nent v. The linear resistivity is nonzero at any finite tem- 
peratures but nonlinear behavior sets in at a crossover 
current with a temperature dependence determined by 
the exponent v. This zero-temperature transition leads 
to slow relaxation dynamics where the correlation length 
diverges as a power law and the relaxation time diverges 
exponentially as the temperature vanishes. 

Simulations of the relaxation dynamics^ found a be- 
havior analogous to relaxation in supercooled liquids with 
a characteristic dynamic crossover temperature rather 
than an equilibrium transition temperature, which is not 
inconsistent with the zero-temperature transition sce- 
nario. On the other hand, a systematic study by MC 
simulations 13 of a sequence of rational values of / con- 
verging to the irrational frustration, using the vortex rep- 
resentation, found two phase transitions at finite temper- 
atures, a vortex-order transition weakly dependent on 
/ and a vortex pinning transition at much lower tem- 
peratures varying with /, which should correspond to 
the resistive transition. These results are in qualitative 
agreement with MC simulations using the phase repre- 
sentation of the same model^ but different ground states 
were found. 

More recently, MC simulations for the the specific 
heat and relaxation dynamics found an intrinsic finite- 
size effect^. The corresponding scaling analysis sug- 
gested a zero-temperature transition with a critical ex- 
ponent v consistent with the value obtained initially 
from current-voltage scaling^. However, a study of 
the low-temperature configurations for frustrations close 
the irrational value by MC simulations in the vortex 
representation^, find two phase transitions consistent 
with earlier worki 3 -. 

On the experimental side, some results on arrays of su- 
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perconducting grains at irrational frustration 2 ^ are con- 
sistent with the scenario of the zero-temperature resis- 
tive transition but on wire networks^, resistivity scaling 
showed evidence of a transition at finite temperature. Re- 
cently, resistivity scaling suggesting a finite temperature 
transition was also observed in arrays of superconducting 
grains^. 

In view of these conflicting results, it seems useful to 
further investigate the current- voltage scaling for the ar- 
ray at irrational frustration by studying both the non- 
linear and linear resistivity with an improved method^ 7 - 
taking into account the long relaxation times. In fact, as 
found recently, current-voltage scaling turned out to be 
quite reliable in determining the phase-coherence transi- 
tion even for a model with quenched disorder, such as the 
three-dimensional XY-spin glass mode l 17 ' 18 . The main 
question is therefore, if the array at irrational frustration 
displays an equilibrium phase-coherence transition at a 
nonzero critical temperature into a state with vanishing 
linear resistivity or its critical temperature vanishes and 
the linear resistivity is finite at nonzero temperatures. 

In this work, we investigate the resistivity scaling of 
Josephson-junction arrays at a irrational frustration / = 
(3 — y/E)/2, a golden irrational, using a driven MC dy- 
namics in the phase representation introduced recentlyi 7 -. 
The results are consistent with a phase-coherence tran- 
sition scenario where the critical temperature vanishes, 
T c = 0. The linear resistivity is finite at nonzero temper- 
atures but nonlinear behavior sets in at a temperature- 
dependent crossover current determined by the thermal 
critical exponent v. The results agree with earlier simula- 
tions using the resistively shunted-junction model for the 
dynamics of the arrayi^. However, with the present MC 
method we are able to reach much lower temperatures 
and current densities, improving the analysis of resistiv- 
ity scaling and the estimate of the critical exponent v. We 
also argue that the finite-temperature transition found 
in resistivity measurements on wire networks^have been 
obtained in a current regime above the crossover current. 

We consider a two-dimensional Josephson-junction 
square array described by the Hamiltonian 

H = -J J2 cos (^ - 0j - Mi) - ~ 6l +^ W 

<ij> i 

The first term gives the Josephson-coupling energy be- 
tween nearest neighbor grains where line integral of the 
vector potential Ay is constrained to Y]^ A.y — 2nf 
around each plaquette. The second term represents the 
effects of an external driving current density J applied 
in the x direction. When J ^ 0, the total energy is un- 
bounded and the system is out of equilibrium. The lower- 
energy minima occur at phase differences 9i — 9i+ x which 
increases with time t, leading to a net phase slippage 
rate proportional to < d(9i — 9i +x )/dt >, corresponding 
to the voltage Vi t i+ X . We take the frustration parameter 
/ equals an irrational number, / = (3— v5) /2, related to 
the Golden Ratio $ = (1 + y/E)/2 as / = 1 - 1/$. In the 
numerical simulations we use periodic (fluctuating twist) 



boundary conditions on lattices of linear sizes L and cor- 
responding rational approximations $ = F n+ i/ F n , where 
F n are Fibonacci numbers (13, 21, 34, 55), with L = F n . 

To study the current-voltage scaling, we use a driven 
MC dynamics method^ 7 -. The time dependence is ob- 
tained by identifying the MC time as the real time t and 
we set the unit of time dt = 1, corresponding to a com- 
plete MC pass through the lattice. Periodic (fluctuat- 
ing twist) boundary conditions are usedi^. This bound- 
ary condition adds new dynamical variables, u a (a — x 
and y), corresponding to a uniform phase twist between 
nearest-neighbor sites along the principal axis directions 
x and y. A MC step consists of an attempt to change the 
local phase 6i and the phase twists u a by fixed amounts, 
using the Metropolis algorithm. If the change in en- 
ergy is AH, the trial move is accepted with probability 
min{l, exp(— AH/kT)}. The external current density J 
in Eq. Q] biases these changes, leading to a net voltage 
(phase slippage rate) across the system given by 

1 d L 

V =Lj t T,( 9 ^-^ + i,j-u x L), (2) 

3=1 

in arbitrary units. The main advantage of this MC 
method compared with the Langevin dynamics used 
earlier— is that in principle much longer time scales can 
be accessed which allows one to obtain reliable data at 
much lower temperatures and current densities. We have 
determined the electric field E = V/L and nonlinear re- 
sistivity p = E/J as a function of the driving current 
density J for different temperatures T and different sys- 
tem sizes. We used typically 2 x 10 5 MC steps to reach the 
nonequilibrium steady state at finite current and equal 
time steps to perform time averages with and additional 
average over 4 — 6 independent runs. 

We have also determined the linear resistivity, pl = 
limj-x) -E/J, from equilibrium MC simulations. As any 
transport coefficient, this quantity can be obtained from 
equilibrium fluctuations and therefore can be calculated 
in absence of an imposing driving current (J = 0). From 
Kubo formula, the linear resistivity (resistance in two 
dimensions) is given in terms of the equilibrium voltage 
autocorrelation as 

PL = hj <ww°)>- ( 3 ) 

Since the total voltage V is related to the phase difference 
across the system A9(t) by V — dA9(t)/dt, we find more 
convenient to determine ph from the long-time equilib- 
rium fluctuations 11 of A9(t) as 

PL= 2^<(A0(*)-A0(O)) 2 ), (4) 

which is valid for sufficiently long times t. To insure that 
only equilibrium fluctuations are considered, the calcu- 
lations were performed in two steps. First, simulations 
using the exchange MC method (parallel tempering) 20 
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FIG. 1: Nonlinear resistivity E/J at different temperatures 
T for system size L = 55. 



were used to obtain equilibrium configurations of the sys- 
tems at different temperatures 21 . This method is known 
to reduce significantly the critical slowing down near the 
transition allowing fully equilibration in finite small sys- 
tem sizes. These configurations were then used as ini- 
tial states for the driven MC dynamics process described 
above, with J = 0, in order to obtain the pl- The ini- 
tial states are similar to the low-temperature states ob- 
tained previousl y 13 ! 16 including thermal excitations. In 
the parallel-tempering method 2 ^, many replicas of the 
system with different temperatures arc simulated simul- 
taneously and the corresponding configurations are al- 
lowed to be exchanged with a probability satisfying de- 
tailed balance. The equilibration time can be measured 
as the average number of MC steps required for each 
replica to travel over the whole temperature range. We 
used typically 4 x 10 6 (parallel tempering) MC steps for 
equilibration which is much larger than the estimated 
equilibration time in systems with up to 100 replicas. 
Subsequent MC simulations for the linear resistivity ob- 
tained from Eq. [4] were performed using 2 x 10 3 time 
averages for 2 x 10 5 MC steps which is much larger than 
the equilibrium relaxation time. 

Fig. la shows the nonlinear resistivity E/J as a func- 
tion of temperature for the largest system size. At small 
current densities J, the nonlinear resistivity E/J tends 
to a constant value, corresponding to the linear resis- 
tivity pl, which decreases rapidly with decreasing tem- 
perature. For increasing J, the resistivity cross over to 
a nonlinear behavior at a characteristic current density 
J n i, which also decreases with decreasing temperature. 
To verify that the nonzero values approached at low cur- 
rents in Fig. 1 correspond indeed to the linear resistivity 
Pl, we show in Fig. 2 the temperature dependence of 
Pl obtained without current bias from Eq.(f5]) for dif- 
ferent system sizes, pl decreases with system size but 
approaches nonzero values for the largest system size. 
These values are in agreement with the corresponding 
values at the lowest current in Fig. 1. Since the be- 



FIG. 2: Temperature dependence of the linear resistivity for 
different system sizes. 



havior of the pl for the largest system size on the log- 
linear plot in Fig. 2 is a straight line, it indicates an 
activated Arrhenius behavior, where the linear resistiv- 
ity decreases exponentially with the inverse of tempera- 
ture with a temperature-independent energy barrier, es- 
timated as Eb ~ 1.07. Such activated behavior suggests 
that the linear resistivity can be very small at low tem- 
peratures but nevertheless remains finite for all tempera- 
tures T > and therefore there is no resistive transition 
at finite temperatures. However, as will be described 
below, the system behaves as if a resistive transition 
occurs at zero temperature, corresponding to a phase- 
coherence transition where the critical temperature van- 
ishes, T c = 0. 

The behavior of the linear resistivity can be related to 
the equilibrium relaxation time for phase fluctuations. 
Since the voltage is the rate of change of the phase, 
a nonzero pl requires measurements over a time scale 
t oc 1/ pl, corresponding to the relaxation time for phase 
fluctuations. Thus, we expect that r should also have 
an activated behavior, increasing exponentially with the 
inverse of temperature. To verify this behavior, we have 
in addition calculated the relaxation time r for different 
temperatures from the autocorrelation function of phase 
fluctuations C{t) as 



1 



C(0) 2 



dtC(t) 



(5) 



using MC simulations with J = 0. The starting con- 
figurations were taken from equilibrium configurations 
obtained 21 with the parallel tempering MC method 20 . 
The results shown on the log-linear plot in Fig. 3 are in- 
deed consistent with an activated behavior of r with an 
energy barrier Ef, = 1.18 in reasonable agreement with 
the value obtained for the linear resistivity in Fig. 2. 

The behavior in Figs. 1, 2 and 3 has the main features 
associated with a phase transition that only occurs at 
zero temperature, T c = 0, similar to the two-dimensional 
gauge glass model of disordered superconductors^^. In 
this case the correlation length £ is finite for T > but it 
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increases with decreasing temperature as £ oc with 
v a critical exponent. The divergent correlation length 
near the transition determines both the linear an nonlin- 
ear resistivity behavior leading to current-voltage scaling 
sufficiently close to the critical temperature and suffi- 
ciently small driving current. To understand in detail 
the behavior of the linear pl and nonlinear resistivity p 
we need a scaling theory for the resistive behavior. If the 
data satisfy such scaling behavior for different driving 
currents and temperatures, the critical temperature and 
critical exponents of the underlying equilibrium transi- 
tion at J = can then be determined from the best data 
collapse. A detailed scaling theory has been described 
in the context of the current-voltage characteristics of 
vortex-glass models^ but the arguments should also ap- 
ply to the present case. The basic assumption is the 
existence of a second order phase transition. Measurable 
quantities should then scale with the diverging correla- 
tion length £ oc |T — T c | _l/ and relaxation time r near the 
critical point. The nonlinear resistivity E/J should then 
satisfy the scaling formi 



T j T = 9±{ jtO' 



(6) 



in two-dimensions, where g± (x) is a scaling function. The 
+ and — signs correspond to T > T c and T < T c , respec- 
tively. If T c 7^ 0, then to satisfy such scaling form, the 
nonlinear resistivity curves on the log-log plot in Fig. I 
should have a positive curvature at small J, with E/J de- 
creasing with decreasing J to a temperature dependent 
value for T > T c while for T < T c , the curvature should 
be negative, with E/J vanishing in the limit J — > 0. The 
data in Fig. 1 do not show a change in curvature even for 
the lowest temperature, already suggesting the possibil- 
ity of a resistive transition at much lower temperatures 
or at T c — 0. However, a full scaling analysis of the data 
is required to show that a transition indeed occur with 
T c = 0. If T c = 0, then the correlation length £ cx T~ v 
and the linear resistivity pl are both finite at T > 0. 
One can then consider the behavior of the dimensionless 
ratio E/J pl which should satisfy the scaling form 



E 

J PL 



(7) 



where g is a scaling function with g(0) = 1. A crossover 
from linear behavior, when g(x) ~ 1, to nonlinear behav- 
ior, when g{x) » 1, occurs when x ~ 1 which leads to a 
characteristic current density at which nonlinear behav- 
ior sets in decreasing with temperatures as a power law, 
J n i oc T/£ oc T 1+u . The scaling form in Eq. ©contains a 
single critical exponent v and does not depend on the par- 
ticular form assumed for the divergence of the relaxation 
time t. However, for sufficiently low temperatures, the 
relaxation process is expected to be thermally activated^ 
with r oc exp(Eb/kT). This corresponds formally to a 
dynamic exponent z — > oo, if power-law behavior is as- 
sumed for the relaxation time x oc From the scal- 
ing form of Eq. © , the linear resistivity should scale as 




FIG. 3: Temperature dependence of the relaxation time r of 
phase fluctuations for system size L = 55. 
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FIG. 4: Scaling plot of the nonlinear resistivity in Fig. 1 for 
v = 1.4. 



Pl oc 1/t and therefore it is also expected to have an 
activated behavior, r oc exp(— Eb/kT). In general, the 
energy barrier Eb also scales with the correlation length 
as Eb oc £V , which leads to a temperature-dependent bar- 
rier Eb oc T~"^ v . A pure Arrhenius behavior corresponds 
to ip = 0. The behavior of the nonlinear and linear resis- 
tivity in Figs 1, 2 and the relaxation time in Fig. 3 are 
quite consistent with these predictions from the scaling 
theory of a zero-temperature transition. 

If there is a zero-temperature transition, as suggested 
by the behaviors in Figs. 1, 2 and 3, then the data for 
the nonlinear resistivity should satisfy the scaling form 
of Eq.([7]), if finite-size effects are negligible, and the best 
data collapse provides an estimate of the critical expo- 
nent v. We expect that finite-size effects are negligible 
for the largest system size L = 55 in Fig. 1 since at 
this length scale the behavior of the linear resistivity is 
roughly independent of the size as can be seen from Fig. 
2. Fig. 4 shows that indeed the data for the largest 
system size satisfy this scaling form with v ~ 1.4 ± 0.2. 

The nonlinear resistivity should also satisfy the ex- 
pected finite-size behavior in smaller system sizes when 
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the correlation length £ approaches the system size L. 
According to finite-size scaling, the scaling function in 
Eq. ([7|), should also depend on the dimensionless ratio 
L /£ and so to account for finite-size effects the nonlinear 
resistivity should satisfy the scaling form 



E 

J PL 



J 



, L X I V T). 



(8) 



The scaling analysis of the whole nonlinear resistivity 
data is rather complicated in this case since the scal- 
ing function depends on two variables. To simplify 
the analysis^ we first estimate the temperature and 
finite-size behavior of the crossover current density J n i 
where nonlinear behavior sets in as the value of J where 
E/JpL = C, a constant. Then, from Eq. ©, the finite- 
size behavior of J n i can be expressed in the scaling form 



JmL^/" = -g(L 1/ "T). 



(9) 



The best data collapse according to the scaling in Eq. © 
provides an alternative estimate of the critical exponent 
v. Fig. 5 shows that indeed the values of J n i for different 
system sizes and temperatures satisfy this scaling form 
with v ~ 1.4, in agreement with the estimate obtained 
for the largest system in Fig. 4 size using Eq. ||7J). 

In addition to the standard finite-size effects, which oc- 
cur when the correlation length is comparable to the sys- 
tem size, already taken into account in the scaling form 
of Eq. {SJ, there are also intrinsic finite-size effects^ 
resulting from the rational approximations used for the 
irrational value of /. Since we use rational approxima- 
tions <f> = F n+ i/F n , where F n are Fibonacci numbers 
(13,21,34,55), with the system size set to L = F n) this 
amounts essentially to have different values of the frus- 
tration, Jl = 1 — 1/$, for different system sizes which 
will only converge to the correct value / = (3 — v / 5)/2 
in the infinite-size limit. We have assumed that such ef- 
fects are negligible in the above scaling analysis but they 
should affect our estimate of the critical exponent v. In 
principle, this intrinsic effect could be taken into account 
within the zero-temperature transition scenario by allow- 
ing for a size-dependent critical temperature T C (L) in the 
scaling analysis 15 . Alternatively, we could regard it as a 
crossover from the critical behavior at the true irrational 
frustration (infinite-size limit) to a phase with an addi- 
tional small frustration Sf = Jl — f which should act as 
a relevant perturbation. In this case, the scaling func- 
tion in Eq. ((7|) should also depend on the dimensionless 
ratio £ 2 5f and again a scaling analysis with more than 
one variable is required. However, our present numeri- 
cal data is not sufficiently accurate to separate this effect 
from standard finite-size effects. 

The present results for the linear and nonlinear resis- 
tivity of the array at irrational frustration obtained by 
the driven MC dynamics agree with earlier simulations of 
the current- voltage scaling using the resistively shunted- 
junction model for the dynamics of the arrayi^, where a 
zero-temperature resistive transition was suggested and 




TL'" 



FIG. 5: Finite-size scaling plot of the crossover current den- 
sity J n i with v = 1.4, for different system sizes L. 



the critical exponent was estimated as v = 0.9(2). Al- 
though the later model is expected to be a more realistic 
description for the dynamics of the array, the value of the 
static critical exponent v should be the same for both 
models. In general, the dynamic exponent z may de- 
pend on the particular dynamics but since the relaxation 
time t is found to diverge exponentially for decreasing 
temperature it corresponds to z — > oo for both dynam- 
ics. The present estimate of v = 1.4(2), however, should 
be more reliable since it considers much lower tempera- 
tures and current densities and larger system size. In- 
terestingly, similar behavior for the resistive transition 
has been found both numerically and experimentally for 
two-dimensional disordered superconductors in a mag- 
netic field described as a gauge-glass modei^i but with 
different value for critical exponent v ~ 2. It should 
be noted however that the actual ground state at irra- 
tional frustration (without disorder) can be quite differ- 
ent, as the self similar structure which has already been 
proposed^ii. As would be expected, the different nature 
of ground state leads to the different values of the critical 
exponent v. 

Although the above scaling analysis is consistent with a 
zero temperature transition, on pure numerical grounds 
the data in Figs 1 and 2 can not complete ruled out 
a vortex-order or a phase-coherence transition at tem- 
peratures much lower than T = 0.15. In fact, phase- 
coherence transitions were found in MC simulations using 
the Coulomb-gas presentation^ at temperatures as low 
as T ~ 0.03 for the sequence of rational approximations 
ft of the irrational / but since they show considerable 
variation with fi, it is not clear if it will remain nonzero 
in the large size limit. However, the lowest temperature 
in Figs 1 and 2 is already much smaller than the apparent 
freezing temperature Tf ~ 0.25 observed in earlier MC 
simulations 8 . Below Tt, a nonzero Edwards- Anderson 
order parameter q(t) —< Si > 2 , was observed, where 
S = (cos 8, sin 8) and the average was taken over the sim- 
ulation times t. Although this could suggest a diverging 
relaxation time r oc J°° q{t)dt near a finite temperature 
T c ~ Tf, such long relaxation time can also result from a 
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zero-temperature transition (T c — 0) as suggested by the 
above scaling analysis since in this case r diverges expo- 
nentially with decreasing temperature, r oc exp(EbfkT), 
as shown in Fig. 3. For low enough temperatures, t will 
eventually be larger than any simulation or experimental 
measuring time scale and an apparent (time dependent) 
freezing transition could occur depending on the partic- 
ular dynamics and system size. 

Some experimental results on arrays of superconduct- 
ing grains at irrational frustration 2 ^ are consistent with 
the scenario of a zero-temperature resistive transition 
since even at the lowest temperatures a zero-resistance 
state was not observed in these experiments. On the 
other hand, current-voltage scaling analysis of experi- 
mental data on wire networks 5 -^ was found to be con- 
sistent with a resistive transition at finite temperature. 
We note, however, that although the equilibrium behav- 
ior of wire networks can be described by the same model 
of Eq. [TJ the nonlinear dynamical behavior may be quite 
different since the nodes of the network are connected 
by continuous superconducting wires, instead of weak 
links, leading to additional larger energy barriers for vor- 
tex motion, not included in the model, and consequently 
larger phase-coherence length £ and relaxation time r 
when compared with weak links 2 ^. In this case, the 
characteristic crossover current to the linear resistivity 
regime at low temperatures due to thermal fluctuations, 
J n i oc kT/£, expected in the zero-temperature transition 
scenario, may only occur at current scales too small to 
be detected experimentally. Thus the resistive behavior 
is observed in a current regime at higher currents where 
it follows the mean-field theory resultS 5 . where a vortex- 
glass transition is possible at finite temperatures. How- 
ever, the zero-temperature resistive could in principle be 
observed in specially prepared wire networks in the weak 
coupling regime where the additional energy barrier for 
vortex motion can be minimized 26 . Other effects, such 
as weak disorder, which is inevitably present in both ex- 
perimental systems, should also be considered. It could 
provide a possible explanation for the finite-temperature 
resistive transition observed recently in arrays of super- 
conducting grains^. 

In conclusion, we have investigated the resistivity scal- 
ing of Josephson-junction arrays at a irrational frustra- 



tion using a driven MC dynamics^!. The results are con- 
sistent with a phase-coherence transition scenario where 
the critical temperature vanishes, T c = 0. The linear 
resistivity is finite at nonzero temperatures but nonlin- 
ear behavior sets in at a crossover current determined 
by the thermal critical exponent v. The results agree 
with earlier simulations using the resistively shunted- 
junction model for the dynamics of the arrayi - and more 
recent MC simulations taking into account the intrinsic 
finite-size effect^. Although we have only studied the 
array at a particular value of irrational frustration, the 
golden mean, we believe that the conclusion of a zero- 
temperature phase-coherence transition should be valid 
for all irrationals but possibly with different values of 
the thermal critical exponent v. The main advantage of 
studying the golden mean value is that it is considered 
the farthest from the low-order rationals and so intrinsic 
finite-effects should be smaller. However, other irrational 
frustrations have also been studied numericall y 1 ^ 23 and 
experimentally 5 . The resistive behavior probes mainly 
the phase-coherence of the system and since we find that 
phase coherence is only attained at zero temperature, we 
can not address directly the question of the existence of 
a vortex-order transition at finite temperatures. In fact, 
vortex order does not require long-range phase coherence. 
Therefore, a vortex-order transition at zero temperature 
or at finite temperature is consistent with the present 
work. However, in view of the results for the supercooled 
relaxation^ 2 - suggesting an analogy to structural glasses 
such transition may be expected at finite temperature 
and in fact is consistent with MC simulations indicating 
a first-order vortex transitio n 13 i 14 i 16 . Thus, the inter- 
esting possibility arises where the array undergoes two 
transitions for decreasing temperature, a finite-resistance 
vortex-order transition at finite temperature and a su- 
perconducting transition only at zero temperature. This 
phase transition scenario and the predicted behavior of 
the linear and nonlinear resistivity provides an interesting 
experimental signature for a Josephson-junction array at 
irrational frustration. 
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